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ABSTRACT 

Statistical regression equations were derived to predict 
future 24-hour changes in intensity of tropical storms and 
typhoons in the western North Pacific. The predictors were 
chosen from 55 parameters available at six-hourly observa- 
tions of tropical storms and typhoons during the period 1960 
-1969. The dependent data were composited into six catego- 
ries: east-west moving storms, recurving storms, and all 

storms within latitude bands 0-9. 9N, 10-19. 9N, 20-29. 9N and 

30-39 .9N. 

The forecast equations were evaluated on a five-year 
(J955-1959) sample of independent data and compared to the 
forecast verification s cheme emp loyed by Fleet Weather Cen — 
ter/Joint Typhoon Warning Center. Two five-predictor equa- 
tions, vjhich require only 12 hours of liistory, can predict 
intensity for the majority of storms within the period July 
-September, and give significantly better results than cur- 
rent intensity forecast methods. 
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INTRODUCTION 



I . 

Intensity (maximum vind) forecasts are one of the more 
difficult aspects of present-day typhoon forecasting. The 
processes governing the intensification of tropical cyclones 
are very complex and not well understood; consequently, 
there is relatively little information available to aid the 
forecaster. Historically, more emphasis has been placed on 
the storm-movement forecast rather than storm intensity, 
since an intensity forecast was of little value to the local 
forecast if the storm-track prediction was inaccurate. Only 
recently, with the increased accuracy of track forecasts, 
has attention turned to intensity. The Joint Typhoon Warn- 
ing Center, Guam (FWC/JTWC), did not attempt to verify in- 
tensity forecasts prior to 1969. The basic intensity fore- 
casting technique used by FWC/JTWC is a linear extrapolation 
of past change in intensity subjectively modified by expected 
conditions along the predicted track (FWC/JTWC 1970). 

A diagnostic approach can be applied to intensity fore- 
casting of tropical cyclones by considering such factors 
as sea-surface temperature, characteristics of the warm air 
mass in which the storm is embedded, and synoptic features 
expected along the predicted storm track. This approach 
mainly indicates whether a tropical cyclone will intensify 
or weaken as it proceeds along the predicted track. 
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Current emphasis is being placed on climatological 
studies of tropical cyclone intensification, and these 
studies are essentially the only quantitative guidance 
available to the intensity forecaster. Gray (1970) has 
provided a general updated climatology and physical theory 
of tropical cyclones in the western North Pacific, A study 
based on 25 years (1945-1969) of data by Brand and Gaya 
(1971) offers geographic and seasonal variations of tropi- 
cal cyclone intensity changes for 12, 24 and 48 hours. 

More recently Liechty (1972) has demonstrated monthly 
variations of intensity of tropical cyclones. 

Although climatological guides are now available to 
aid the forecaster in predicting the future intensity of 
tropical cyclones, the purpose of this study is to provide 
the forecaster with an additional tool to forecast the 
future 24-hour intensity of tropical cyclones in the western 
Pacific, This tool will be in the form of prediction equa- 
tions, derived using statistical regression techniques, 
which will quantitatively predict the future 24-hour change 
in intensity of tropical cyclones. 
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II. DATA 



A. DATA SOURCE 

The data used in this study were obtained from the Naval 
Environmental Prediction Research Facility, Monterey, Calif- 
ornia. The data source consisted of a history file of six 

1 2 

hourly information on tropical storms and typhoons (here- 
after referred to as storms) which occurred in the western 
North Pacific Ocean during the period 1945 through 1969 and 
are based on post-season analysis. These data were compiled 
by the National Climatic Center, Asheville, North Carolina, 
for the Navy Weather Research Facility, Norfolk, Virginia, 
and have served as the basis f o.r several studies on intensity 
of tropical storms and typhoons (Brand and Gaya 1971). More 
recently, the same data were used in a study of intensity 
changes of tropical cyclones for the period 1960-1969 
(Liechty 1972). 

B. DEPENDENT DATA SAMPLE CRITERIA 

Data were chosen for the dependent sample from the years 
1960 through 1969, since it was felt that the last few years 
of wind data would be the most complete and accurate consi- 
dering the increased use of Doppler radar. Following Liechty 
(1972), it was decided to eliminate storms of the following 
categories : 

^Tropical cyclonic circulation which attains tropical storm 
intensity (34-63kt) at one time in the life of the storm. 

2 

Tropical cyclonic circulation v;hich attains typhoon inten- 
sity (^ 64kt) at one time in the life of the storm. 
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!• Storms with an initial northeasterly heading which 
continued to move to the northeast. 

2. Storms which had a genesis west of 125E. 

3. Looping storms, i.e. storms which formed a complete 
loop in the track. 

4. Storms with the first observation reporting a maxi- 
mum intensity greater than 65 knots. 

In addition, each tropical cyclone was required to have 

t 

a 48-hour history. Due to the strong effects of terrain on 
the behavior of tropical cyclones, as suggested by Brand 
(1972), only those storms with continuous track segments 
over open ocean were considered. V7hen a storm track reached 
land, the data were discontinued at the preceding six-hourly 
observation time. 

Considering the large variability in the storm tracks, 
the data were first classified into two categories: 

1. East-west moving storms, i.e. storms moving pre- 
dominately in a westerly direction and not looping or 
recurving . 

2. Recurving storms, i.e. storms which changed direc“ 
tion of movement in a clockwise direction from a westerly 
to an easterly component.- 

The data were further classified into monthly periods, 
since tropical cyclones exhibit significant monthly and 
seasonal variations in intensity (maximum surface wind) (see 
Liechty 1972), Table I is a summary of the dependent data 
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TABLE I 



Dependent data summary with means and standard 
deviations of future 24-hour change in intensity. 





(a) East 


-VJest Moving Storms 




Period 


Number of 
Storms 


Number of 
Six-hourly 
Observations 


Mean 


Std . 
Dev . 


July 


18 


250 


1.8 


23.0 


August 


22 


179 


9 . 8 


26 . 2 


September 


25 


319 


4.6 


21 . 3 


October 


8 


75 


9.1 


22.6 


November 


8 


102 


6 . 9 


19 . 3 


De c- Jun 


17 


149 


5.3 


23 . 1 


Jul-Aug-Sep 


65 


748 


4.9 


23 . 3 




(b) Recurving Storms 






July 


2 


41 


-1.5 


11.5 


Augus t 


12 


250 


0.9 


19 . 4 


September 


13 


321 


3.0 


20 . 1 


October 


■ 20 


380 


1.5 


22.5 


November 


9 


206 


1.5 


22.3 


Dec- Jun 


12 


246 


-0.4 


23.6 


Jul-Aug- Sep 


27 


612 


1.8 


19.4 
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sample grouped by months in the east-west and recurving 
categories • 

The data included storms which both intensified and 
weakened and many which went through a complete cycle of 
intensification and weakening. Since all storms were only 
considered over open ocean, all of the track over land was 
eliminated, thus biasing the sample somewhat in favor of 
storms which intensified. The magnitude of the means of 
the future 24-hour change in intensity was in general lar- 
ger in the east-west moving storms than in the recurving 
storms, and all mean values were positive in the east-west 
storms. The largest mean values occurred in August and 
September storms. The standard deviations were slightly 
larger in east-west moving storms; hov7evcr, the standard 
deviations of both east-west and recurving storms are very 
similar in magnitude. The largest standard deviation of 
future 24-hour change in intensity occurred in August east- 
west moving storms. These means and standard deviations 
cannot be compared to the work done by Liechty (1972) since 
he was concerned with actual intensities rather than 
changes. Since all observations were over open ocean, the 
weakening of many east-west moving storms was not seen in 
the data, due to discontinuing data 24-hours prior to land- 
fall, However, the entire track of most recurving storms 
was observed, thus accounting for a much larger mean value 
of future 24-hour change in intensity of east-west moving 
storms. 
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In the western Pacific, the frequency of occurrence of 
typhoons is correlated with potential instability and water 
temperature. Temperature is thought to be an important fac 
tor in tropical cyclone intensity, and intense tropical 
cyclones are confined to areas with warm surface waters 
(Palmen and Newton, 1969) • In view of the fact that temper 
ature information was not available in this data, the data 
were also stratified into four latitude bands from 0-9.9N, 
10-19*9N', 20-29 * 9N and 30-39. 9N. Stratification of data in 

this manner, to a certain extent, took into account the 
variation of sea-surface temperature with latitude in the 
Pacific. The data stratified by latitude bands were classi- 
fied into monthly periods. Table II is a summary of depen- 
dent data stratified into latitude bands. 

When examined in latitude bands, the mean values of 
future 24-hour change in intensity were positive in the 
lov7er latitudes, 0-20N. At higher latitudes, 20-40N, the 
mean value of future 24-hour change in intensity decreased 
to a small positive value in August and negative values in 
all other months. The largest mean value was observed in 
intensifying August storms, in the 10-19. 9N latitude band, 
and in the V7eakening Nov-Dec storms within 20-29. 9N, As 
storms moved northward out of warm southern waters, they 
would be expected to weaken, and this was particularly 
evident in the winter and spring storms (Nov-Jun) moving 
from 10-19. 9N to 20-29. 9N. 
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TABLE II 



Dependent data summary in latitude bands with means 
and standard deviations of future 24-hour change in 
intensity . 

(a) 0-9. 9N 



Period 



Number of 
Six-hourly 
Observations 



Mean 



Standard 
Devia t ion 



Jan-Dec 


190 




11.2 


19 . 5 




(b) 


10-19 . 9N 






Dec- Jun 


251 




3.3 


23.0 


July 


136 




10.6 


19 . 9 


Aug us t 


122 




15 . 2 


22 . 9 


September 


288 




12 . 7 


21 . 0 


October 


227 




13.3 


20.1 


November 


197 




6.7 


18.1 


Jul-Aug-Sep 


546 




12 . 7 


21 . 2 




(c) 


20-29 . 9N 






Nov - J un 


90 




-22,1 


19 . 0 


J uly 


132 




-5.9 


20 . 7 


Au gust 


211 




2.3 


23 . 2 


September 


330 




-3.5 


17.0 


October 


217 




-7.1 


20.1 


Jul-Aug-Sep 


673 




-2 . 2 


20.1 




(d) 


30-39 . 9N 






Jan-Dec 


127 




-12.2 


12 . 3 
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The standard deviation of future 24-hour change in 
intensity was generally smaller when categorized into lati- 
tude bands, thus indicating the importance of subtracting 
out different means in different latitudes. However, the 
standard deviations were still rather uniform. 

C. INDEPENDENT DATA SAMPLE CRITERIA 

The independent data sample used for equation verifica- 
tion were chosen from the years 1955 through 1959. Storms 
in the months of July, August and September were chosen for 
verification. Once again, only observations of storms over 
open ocean were used. The data were not divided into east- 
west or recurving categories, because it would be impossible 
to know in the early stages if a storm was going to recurve. 
Storms of the following categories were eliminated: 

1. Northeasterly moving storms v/hich had an initial 
northeasterly heading. 

2. Storms which had a genesis west of 125E, 

3. Storms with the first observation reporting a 
maximum intensity greater than 65 knots. 

4. Looping storms . 

Each tropical cyclone was required to have a 48-hour 
history, and data were discarded 24 hours prior to landfall, 
to allow 24-hour forecasts only over open ocean. 

The data were considered in monthly periods, and were 
also stratified into four latitude bands. Table III is a 
summary of the independent test data. The mean values of 
the future 24-hour change in intensity of the independent 
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TABLE III 



Independent data summary by month with means and 
standard deviations of future 24- hour change in 
intensity. 

(a) Monthly Category 
Number of 



Period 


Number of 
Storms 


S ix-hour ly 
Observations 


Mean 


Standard 
Deviat ion 


July 


11 




174 


1.9 


23.1 


Au gus t 


12 




158 


6.1 


19 . 8 


September 


16 




253 


8.3 


23.4 






(b) 


Latitude Band 


0-9 . 9N* 








(c) 


Latitude Band 


10-19 . 9N 




July 






69 


18.6 


15 . 0 


Augus t 






75 


17.5 


19 . 3 


September 






140 


15 . 1 


20 . 6 






(d) 


Latitude Band 


20-29 . 9N 




July 






92 


-9 . 3 


21.3 


Augus t 






70 


-2.4 


18 . 5 


September 






112 


-0.3 


24 . 1 






(e) 


Latitude Band 


30-39 . 9N* 





^ Indicates that this Latitude Band was not evaluated due 
to small number of observations. There was 1 observation 
within 0-9 . 9N and 26 observations within 30-39. 9N. 
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sample were negative at higher latitudes, just as in the 
dependent sample data. Likewise, the standard deviations 
were in general smaller when stratified into latitude 
bands, just as in the dependent sample. The largest mean 
values again occurred in August and September storms. The 
standard deviation of the future 24-hour change in inten- 
sity of the independent sample is similar to the standard 
deviation of the dependent sample. 

D. DESIGNATION OF DATA 

The variables available for this study are listed in 
Table IV. Variables Z through were obtained as 

1 1 u 

original data from the Environmental Prediction Research 
Facility. These parameters were originally derived for use 
as analog data elements in a study of typhoon movement 
prediction through analogs by Jarrell and Somervell (1970). 
Variables . Z^^ through Z ^2 derived to have available 

past changes of intensity, latitude and longitude. Variables 
Z ^2 through Z^g were computed as "perfect forecast" para- 
meters. Several -of the parameters were combined or manipu- 
lated to describe large scale synoptic factors related to 
storm location and intensity. Z^^ and Z^^^ were used to 
generate Z^^ which is a gradient value of 700mb ridge height 
difference divided by latitude separation of the 700mb ridge 
and the storm. Z^^ and Z^^ were used to generate Z^^ which 
is a gradient measure of 700mb trough height difference 
divided by the longitude separation of the 700mb trough and 
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the storm. A number of predictors were formed by ratios 
and differences of parameters which Liechty (1972) related 
to storm intensity. (intensity minus minimum sea-level 

pressure) is such a relation of intensity and sea-level 
pressure, (maximum intensity divided by storm size) 

is a measure of the relative strength and size of the storm. 
Storm size is defined as the average radius, in degrees 
latitude, of the outer closed isobar. Variables Z^^ and 
Z 48 used to generate Z^^ which is a relative measure 

of the past 12-hour changes of intensity and size. To avoid 
dividing by zero, if the past 12-hour change in size were 
zero, the constant 10.0 was added to and likewise added 

to Z^g. Parameter Z^g (square of maximum intensity) was 
generated to provide a second order variable, and 
(square of intensity divided by size) was designed as a 
relative measure of intensity with respect to size. Although 
sea-level pressure and intensity are highly correlated, Z^^ 
(intensity divided by sea-level pressure) was created as a 
relative measure of intensity and sea-level pressure. A 
measure of the speed of movement related to intensity of 
the storm was formed as Z^^ (24-hour speed divided by 
intensity). Z^g (700mb rainimum height minus sea-level 
pressure) was a measure of thickness of the storm. Z^^ 

(past 24-hour change in intensity divided by size) was 
designed as a relative measure of intensification and size 
of a storm. 



19 



TABLE IV 



Sea-level (and 700MB where specified) predictors 
available at each six-hourly observation. 



Z . 

1 



Variable 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 



Latitude-degrees and tenths 
Longitude-degrees and tenths 
Past 12-hour direction 
Past 12-hour speed, in knots 
Past 24-hour direction 
Past 24-hour speed 

Size-the average radius of the outer closed 
isobar, in degrees latitude 

Past 12-hour change in size, in degrees latitude 
Minimum sea-level pressure, in mbs 
Past 12-hour change in sea-level pressure, 
in mb s 

Maximum intensity-maximum surface wind speed, in 
kno t s 

Minimum 700mb height in tens of meters 
700mb ridge la t i t ud e -no r t h of storm, in whole 
degrees 

700mb ridge height-at ridge line north 
in tens of meters 

700mb trough longitude-at 35 N-nearest 
west of storm, in whole degrees 
700mb trough height-at intersection of 
line at 35°N in tens of meters 
Past 12-hour change in intensity 
Past 24-hour change in intensity 
Past 12-hour latitude change 
Past 12-hour longitude change 
Past 24-hour latitude change 
Past 24-hour longitude change 
Future 12-hour latitude change 
Future 12-hour longitude change 
Future 12-hour 700mb ridge latitude change 
Future 12-hour 700mb ridge height change 
Future 12-hour 700mb trough longitude change 
Future 12-hour 700mb trough height change 
Future 12-hour max intensity change 
Future 24-hour latitude change 
Future 24-hour longitude change 
Future 24-hour max intensity change 
Future 24-hour 700mb ridge latitude change 
Future 24-hour 700mb ridge height change 
Future 24-hour 700mb trough longitude change 
Future 24-hour 700mb trough height change 
Future 12-hour sea-level pressure change 
Future 24-hour sea-level pressure change 



of storm 
t rough 
trough 



20 



z , 

1 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 



TABLE IV 
(Continued) 



Variable 



700mb ridge latitude minus storm latitude 
^^ 13 " ^ 1 ^ 

Max wind speed minus minimum SLP ” ^ 9 ^ 

700mb ridge height minus min 700mb height 
^^14 " ^ 12 ^ 

700mb ridge HT DIFF/LAT SEP of 700mb ridge 
and storm 

700mb trough height minus min 700mb height 
<^16 - 

700mb trough long, minus storm long. 

Trough HT DIFF/trough minus storm long. DIFF 

Max intensity/size ^^xi^^7^ 

Past 12-hour change in size plus 10 (Z^+10.0) 

o 

Past 12-hour change in max intensity plus 10 

2 



(Z^^+10.0) 



Square of max intensity 
Square of max in t ens i t y / s i ze (Z, /Z-) 

Max int ens i ty /SLP (Z^^/Z^) 

Past 24-110111* speed/max intensity 



700mb min height minus SLP (Z 



12 






Past 24-hour change in max int/size (Z_o/Z_) 

1 o / 

Past 12-hour change in max int 4-10.0/Past 
12 -hour change in size + 10.0 (Z^g/Z^^) 
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In this study, emphasis was placed on generating an 
equation to predict ^^2’ 24-hour future change in maxi- 

mum intensity. However, equations to predict actual 24- 
hour future values of intensity were also derived for 
several categories and are presented in Appendix B. 

Results are mainly presented for 24-hour future inten- 
sity changes, rather than actual intensity values. There 
are several papers concerning spatial and seasonal varia- 
tions in 24-hour intensity changes which are available to 
the forecaster, namely Riehl (1971) and Brand and Gaya 
(1971). These papers can be consulted in conjunction with 
the use of the equations from this study. 
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III. PROCEDURES AND METHOD OF ANALYSIS 



A. THE REGRESSION METHOD 

Many problems in research require extensive analysis of 
large amounts of data. The Biomedical Computer Programs 
(Dixon 1966) v;ere developed to provide programs for these 
commonly required tasks of data processing and statistical 
analysis. The program used in this study was the stepwise 
regression analysis program BIMED 02R. 

BIMED 02R computes, in a stepwise manner, a sequence of 
multiple linear regression equations. One variable is 
added to the regression equation at each step. The variable 
added is the one v/hich achieves the greatest reduction in 
the error sum of squares. This variable is also the one 
which has the highest partial correlation with the dependent 
variable at the particular step in the analysis of variance. 
Equivalently, it is the variable which would have the high- 
est F-value if it were added. The F-value or F-statistic 
upon entry, F, is expressed at step k as (Dixon 1966) 

„ Z (C.E.V.k) - Z (C.E.V.k-1) 

(l,n-k-l) % (i.E.V.k) 

where 

%(C.E.V.k) is the percent cumulative explained variance 
at step k 

% (C . E . V . k-l) is the percent cumulative explained 
variance at step k-1 
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%(U,E,V,k) is the percent unexplained variance remaining 
at step k 

The F“Value can be used to determine if a variable is 
statistically significant at a particular level* For 
example, a test at the 5% significance level whether the 
true partial regression coefficients are zero yields a 
critical F-value for particular degrees of freedom. As 
long as the F-value upon entry is larger than the critical 
F-value, it can reasonably be concluded that the variance 
accounted for by regression is greater than could be 
expected if all the true partial regression coefficients 
were zero (Crow, Davis, Maxfield 1955). 

B. RELATED STATISTICAL PARAMETERS 

Several related statistical parameters are available 
as outputs from the BIHED 02R program. Among these para- 
meters are: 

1. multiple correlation coefficient, R 

2. standard error of estimate, S.E. 

3 . mean value , I 

4. standard deviation, O 

5 • F-value 

R is the simple correlation coefficient if only one, 
predictor is added, and R becomes the multiple correlation 
coefficient after k predictors are added. The F-value 
becomes the overall F-value for the equation after k 
predictors are added. The overall F-value of the equation 
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can be compared to a critical F-value, for five predictors 
and a given number of independent samples, for a particular 
significance level to determine if the equation is 
statistically significant. 

C. DERIVATION OF EQUATIONS 

The BIMED 02R program was used to generate equations 

to predict future 24-hour change in intensity 2^2* 

predictors were selected from parameters through Z ^2 

and Z^g through Z^^ in Table IV. Variables Z^^ through 

Z«o were not allowed as predictors since they were future 
J o 

parameters and would perhaps not be available to the fore- 
caster. Thus the prediction equation would not contain 
any parameters other than those available from history or 
from a current aircraft observation. Since variables ’^22 

through Z« were not allowed to enter as predictors in the 
i o 

regression equation, it was not determined if any of these 
parameters were of statistical significance and v;ould 
indeed enter as predictors. 

This study employed two statistical models to generate 
future 24-hour change in intensity I. The first model 
was of the form: 

I=C^+CiN^+C2N2+C3N3+C^N^+C3N5+C^N^+C7N7+C8N8+CgNg+CioNio 

where through and N^ through N^^ were determined 

by the stepwise, least squares technique. through 

are constants and N^ through N^^ are predictors. The units 



25 



of I are knots per day. The second model was of the same 
form as the first, but was limited to five predictors. 
Although some parameters that were still of statistical 
importance were lost in using only five predictors, the 
loss of useful information was very little. The differ- 
ence in results when using five predictors vice ten pre- 
dictors will be discussed later. 

All predictors composing each five-predictor equation 
were tested using the F-value upon entry of the predictor 
against a critical F-value (Crow, Davis and Maxfield 1955) 
and were found to be statistically significant at the 5% 
level. In addition, the overall F-value of each equation 
was compared to a critical F-value for a five-predictor 
equation, and each equation was found to be significant at 
the 0,5% level . 

Regression equations to predict Z ^2 (future 24-hour 
change in intensity) using five predictors were first de- 
rived from the east-west dependent sample data in the 
following seven categories: (a) July (b) August (c) Septem- 

ber (d) October (e) November (f) December through June and 
(g) July through September, A summary of the east-west 
data equations is found in Appendix A, Table X, In the 
east-west equations, variable (square of intensity 

divided by size) entered in five of the seven equations. 
Variables Z^^ (past 12-hour change in SLP) and Z^^^ (past 
12-hour change in intensity) each entered in four of the 
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seven categories. Variables (current latitude) and 

(current longitude) each entered in three of the seven 
categories. 

Regression equations were next derived from the 
recurving dependent data in the same monthly categories 
as in the east-west dependent sample and are in Appendix A, 
Table XI. In the recurving data equations, variable 
(past 12-hour change in SLP) entered in each equation. 

Zj^ (current latitude) entered in six of the seven equations 
Variable (past 12-hour change in intensity) entered in 

four equations . 

If there should be a bias in a current storm and a 
prediction equation could be derived which would detect 
this bias and use it in a new equation as a predictor, the 
results of the new prediction equation could be improved. 

A scheme was tried whereby a ten-predictor equation was 
generated to predict Z (past 24-hour change in intensity) 
The predicted past 24-hour change in intensity derived from 
this equation was then subtracted from the actual past 
24-hour change in intensity, thus generating an error 
parameter to represent a possible bias in the storm. The 
error was then made available as a variable for a ten- 
predictor equation to predict ^^2 (future 24-hour change 
in intensity). This scheme was tried for two cases, Septem 
ber east-west data and August east-west data. In the 
September east-west case, the error parameter did not enter 
into the nev; ten-predictor equation to predict ^ 
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stepwise, least squares technique. In the August east-west 
case, the error parameter did enter as a predictor, but the 
overall results of this new equation were not improved. 

The technique of generating an error to act as a measure of 
bias, and also show how well an equation was performing in 
the past to aid in generating a better equation for future 
predictions, did not prove to be useful in these two cases 
and was not tried further. 

Regression equations to predict (24-hour future 

change in intensity) using five predictors were next derived 
from all data in the latitude band 10-19. 9N for the same 
seven categories as in the east-west and recurving catego- 
ries, and are located in Appendix A, Table XII. In these 
equations, variable (past 12- hour change in SLP) entered 

in six of the seven equations. Variable Z^^ (square of 
current intensity) entered in four of the seven equations. 
Variables Z^ (current latitude), (past 12-hour change 

in intensity), and Z^^Csquare of intensity divided by size) 
each entered in three equations. 

Regression equations to predict Z ^2 (future 24-hour 
change in intensity) using five predictors were next derived 
from all data in the latitude band 20-29. 9N for the follow- 
ing six categories: (a) July, (b) August, (c) September, (d) 

October, (e) July through September and (f) November through 
June. November storms could not be examined separately due 
to the small number of observations in November. 
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The equations to predict Z ^2 20-29. 9N are presented 

in Appendix A, Table XIII. Variable (past 12-hour change 

in SLP) entered in five of the six equations. Variables 

(current latitude), Z^^ (square of current intensity), and 

Z^Q (square of intensity divided by size) each entered in 

three of the six equations. It should be noted that latitude 

entered as a predictor even after separation by latitude, 

suggesting that perhaps a finer degree of separation could 
0 

have been used. 

There was not enough data in each monthly category for 
the latitude bands 0-9 . 9N and 30-39. 9N to justify generating 
a prediction equation for each monthly period. Appendix A, 
Table XIV is a summary of the equations to predict Z^^ 

(future 24-hour change in intensity) using five predictors. 
These equations were derived from all data within 0-9 . 9N 
and 30-39. 9N. 

D. EVALUATION OF EQUATIONS 

1 . Dependent Data Results 

The multiple correlation coefficient, (R) , and the 
standard error of the estimate, (S.E.), of each five predic- 
tor equation were compared to R and S.E. of persistence 
(past 24-hour change in intensity) when used for predicting 
the future 24-hour changes in intensity. A 24-hour forecast 
was made at each 6-hourly observation time using both the 
five-predictor equation and persistence (linear extrapola- 
tion). These forecasts were compared to actual future 
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24-hour changes in intensity, A tabulation of R and S,E. 
for both equation and persistence for each dependent -data 
category is presented in Table V, In each category, the 
correlation coefficient for predicting future 24-hour 
change in intensity was higher for the five-predictor equa- 
tion than for persistence , and correspondingly, the 
standard error of the estimate was lower for the equation 
than for persistence. 

In general, when using equations to predict future 
24-hour changes in intensity, the larger the sample size, 
the smaller the correlation coefficient. The most difficult 
period to predict, i.e., the period with the smaller corre- 
lation coefficient, was the sample from the combined months 
of July, August and September. Using persistence to predict 
future 24-hour changes in intensity, a smaller correlation 
coefficient did not always correspond to a larger sample 
size. Overall, the prediction equation still gave better 
results than persistence, in each category. 

2 • Independent Data Results 

The effectiveness of the statistical regression 
equations in predicting future 24-hour change in intensity 
was first examined as in the dependent data by looking at 
the correlation coefficient (R) and standard error of the 
estimate (S.E.). The predicted change in intensity was 
compared to actual future 24-hour change in intensity. R 
and S.E. for predicting future 24-hour change in intensity 
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TABLE V 



Correlation coefficient (R) and standard error of estimate 



(S.E.) of five predictor equation and persistence for pre- 
dicting future 24-hour change in intensity of dependent 
sample. 



Equation 




Sample 


a. East-West Data 
Equation 


Persistence 


Category 




Size 


R 


S.E. 


R 


S.E. 


JULY 




250 


. 73 


15.9 


.39 


21.3 


AUG 




179 


. 65 


20.1 


. 06 


26.2 


SEPT 




319 


.74 


14.5 


.53 


18 . 1 


OCT 




75 


. 83 


13 . 2 


.26 


22 . 0 


NOV 




102 


.82 


11 . 4 


.53 


16 . 4 


DEC-JUN 




149 


. 68 


17.3 


.29 


22 . 2 


JUL-AUG- 


SEPT 


748 


.66 


17.5 


.37 


21.7 


JULY 




41 


b. Recurving Data 

.94 4.1 


.84 


6 . 3 


AUG 




250 


. 70 


13.9 


.32 


18 . 4 


SEP 




321 


. 71 


14.2 


.45 


18 . 0 


OCT 




380 


. 76 


14.8 


.40 


20 . 6 


NOV 




206 


.80 


13.6 


.48 


19.6 


DEC-JUN 




246 


. 75 


15.7 


.37 


22 . 0 


JUL-AUG- 


SEPT 


612 


.69 


14 . 1 


.42 


17 . 7 


JAN-DEC 




190 


c. 0-9 
.74 


. 9N 

13 . 2 


.42 


17 . 8 


JULY 




136 


d. 10- 
. 78 


19 . 9N 

12 . 7 


. 25 


19 . 4 


AUG 




122 


.78 


14 . 6 


. 33 


21.7 


SEPT 




288 


.64 


16 . 3 


.44 


19 . 0 


OCT 




227 


.68 


14 . 8 


. 31 


19 . 1 


NOV 




197 


. 80 


11 . 1 


.28 


17.5 


DEC-JUN 




251 


. 67 


17 . 3 


.29 


22 . 1 


JUL-AUG- 


SEPT 


546 


. 6 3 


16 . 5 


.37 


19 . 8 


JUL 




132 


e . 20- 

. 69 


2 9 . 9N 

15 . 3 


.37 


19 . 3 


AUG 




211 


. 64 


18.0 


. 05 


23.2 


SEPT 




330 


. 63 


13.3 


.39 


15 . 7 


OCT 




217 


. 80 


12 . 1 


. 26 


19 . 5 


NOV-JUN 




90 


. 83 


10 . 8 


.27 


18 . 4 


JUL-AUG- 


SEPT 


673 


.59 


16 . 2 


.28 


19 '. 3 


JAN-DEC 




127 


f . 30- 

.58 


39. 9N 

10 . 2 


.04 


12.4 
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were compared with corresponding R and S.E. values for a 
forecast using only the past 24-hour change in intensity 
(persistence). These values are presented in Table VI for 
independent data. Just as in the dependent data sample, 
the correlation coefficient of the equation was higher in 
each category than that of persistence when predicting 
future 24-hour change in intensity. In general, equations 
from combined July, August, and September data gave higher 

r 

correlation coefficients than equations derived from one 
month. The most difficult month to forecast was September, 
and the best results were seen in August when examining 
correlation coefficients, September was also the most 
difficult month when using only persistence, and July and 
August were about equal in forecasting with persistence, 

A second, perhaps more timely method of evaluating 
the performance of each equation is the test recently used 
by the Fleet Weather Cen t r al / Jo in 1 Typhoon VJarning Center 
at Guam, Marianas Islands. When considering the standard 
error of estimate (S.E.) and correlation coefficient (R) , 
equal weight is given to the error in the future 24-hour 
change of intensity of a storm, regardless of the actual 
intensity of the storm. However, a 20 knot error in a very 
intense typhoon is not as important as the same error in a 
weak tropical storm. To overcome this deficiency, FWC/JTWC 
suggested describing errors as a fraction of the observed 
wind : 

Error = Forecast Intensity - Observed Intensity 

Observed Intensity 
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Some acceptability criteria (from the viewpoint of adequacy 
for disaster control planning) were established by FWC/JTWC 
for a 24-hour forecast as follows: 



Accurate to within 

measurement error Error ^ 10% 

Adequate Error £ 20% 

Useful Error £ 30% 

Inadequate Error > 30% 

The criteria became less stringent for 48-or 72-hour 
forecasts (FWC/JTWC 1970). The distributions of FWC/JTWC 
1970 and 1971 24-hour intensity forecasts, according to 
acceptability criteria, were as follows (FWC/JTWC 1970 and 
EG PACOM 1972) : 



Accurate to within 


1970 


1971 


measurement error 


31% 


25% 


Ade qua t e 


54% 


55% 


Useful 


7 0% 


7 2% 


Inadequate /Misleading 


30% 


2 8% 



These results provide a guideline for evaluating 
the statistical regression equations using the independent 
sample from 1955-59. For example, a typical error of 10- 
15 knots in 24 hours could represent an accurate forecast 
for an intense typhoon or an inadequate forecast for a weak 
tropical cyclone. 

a. Five-predictor equations derived from east-west 
and recurving data 

A 24-hour forecast was made at each 6-hourly 
observation for data grouped by monthly periods. Since it 
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TABLE VI 



Correlation coefficient (R) and standard error of estimate 
(S.E.) of persistence and five predictor equation for pre- 
dicting future 24~hour change in intensity of independent 
sample . 



Equation 


Sample 

Size 


a • 


July Data 
Equat ion 

R • S • E • 


Persistence 

R t S • E • 


July E-W 


174 




.73 


16 . 0 


. 43 


21.0 


July Rec. 


174 




.53 


19 . 6 


. 43 


21.0 


July 10-19. 9N 


69 




.40 


13 . 8 


. 26 


14.6 


July 20-29. 9N 


92 




.59 


17 . 3 


.40 


19 . 7 


July- Au g- Sep : 
E-W 


174 




.74 


15 . 7 


.43 


21.0 


Rec . 


174 




. 74 


15 . 6 


.43 


21.0 


10-19 . 9N 


69 




.54 


12 . 7 


. 26 


14.6 


20-29 . 9N 


92 




.57 


17 . 6 


.40 


19 . 7 


Aug E-W 


158 


b . 


Au gus t 
. 75 


Data 
14 . 1 


. 42 


19 . 6 


Aug Rec. 


158 




. 79 


13.1 


.42 


19 . 6 


Aug 10-19. 9N 


75 




.52 


16 . 5 


.10 


19 . 3 


Aug 20-29. 9N 


70 




. 85 


9 . 7 


.51 


16 . 1 


Jul-Au g-Sep : 
E-W 


158 




. 82 


12 . 4 


. 42 


19 . 6 


Rec . 


158 




.76 


14 . 0 


. 42 


19 . 6 


10-19 . 9N 


75 




. 62 


15 . 2 


. 10 


19 . 3 


20-29 . 9N 


70 




. 80 


11.1 


.51 


16.1 


Sep E-W 


253 


c . 


September Data 
.56 19.4 


. 17 


23 . 1 


Sep Rec. 


253 




.58 


19 . 0 


. 17 


23 . 1 


Sep 10-19. 9N 


140 




.44 


18.6 


. 19 


20.3 


Sep 20-29 . 9N 


112 




. 68 


17 . 7 


.06 


24.2 


Jul-Aug-Sep : 
E-W 


253 




.66 


17 . 7 


.17 


23.1 


Rec . 


253 




. 63 


18 . 3 


.17 


23 . 1 


10-19 . 9N 


140 




. 51 


17 . 8 


.19 


20 . 3 


20-29 . 9N 


112 




. 83 


13.3 


.06 


24 . 2 
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was unknown whether a storm would be an east-west moving 
storm or a recurving storm, a 24-hour forecast was made 
using equations derived from each category of dependent 
data. Both July east-west and July recurving equations were 
used on July data. In addition, equations derived from the 
combined dependent sample months of July, August and Septem- 
ber in both east-west and recurving categories were used 
for each monthly period , The future 24-hour change in in- 
tensity 'was obtained either from the appropriate equation 
in Appendix A, Table X or XI, or using persistence (past 
24-hour change in intensity) and added to present intensity 
to give the 24-hour forecast intensity. 

Parts a, b and c of Table VII are distributions 
of 24-hour forecasts of intensity using statistical regres- 
sion equations. These are grouped by monthly periods from 
July through September. Part d. of Table VII, is the dis- 
tribution of 24 -hour forecasts using persistence. The fore- 
casts made in Table VII, part d, v?ere not subjectively 
modified in any manner, yet the average distribution by 
persistence corresponds very closely to the distribution 
of the forecasts by FWC/JTWC for 1970 and 1971 as shown 
previously. If samples are similar, this suggests the 
official forecasts in these years were nearly equivalent 
to persistence. In the following comparisons, forecast 
results better than a persistence forecast may then be 
considered better than official forecasts. 
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TABLE VII 



Acceptability distribution of 24-hour intensity forecasts 
for independent data in percent of saraple size, 

a. July Data - Sample Size 174 



Equation : 


July 




J uly 


Jul-Aug-Sep 


Jul-Aug-Sep 




Eas t-Wes t 


Recurving 


East-West 


Recurving 


Accurate 


44 




42 


49 


47 


Adequate 


71 




59 


71 


77 


Useful 


83 




76 


86 


84 


Inadequate 


17 




24 


14 


16 




b. 


August Data - 


Sample Size 


158 


Equa t ion : 


August 




Augus t 


Jul-Aug-Sep 


Jul -Aug-S ep 




East-West 


Recurving 


East -Wes t 


■ Recurving 


Accurate 


40 




52 


53 


54 


Adequate 


79 




82 


82 


83 


Useful 


87 




89 


91 


91 


Inadequate 


13 




11 . 


9 


9 




c . 


September Data 


. - Sample Size 253 


Equation : 


September 


September 


Jul-Aug-Sep 


Jul-Aug-Sep 




East-West 


Recurving 


East -Wes t 


Re cu rv ing 


Accurate 


44 




38 


43 


42 


Adequate 


67 




70 


71 


71 


Useful 


90 




88 


86 


88 


Inadequate 


10 




12 


14 


12 




d . 


24 


-Hour forecast using linear 






extrapolat ion 


(per sis t ence ) 


Period: 


July 




August 


September 


Av e r a g e 


Accurate 


25 




38 


30 


31.1 


Adequate 


50 




66 


51 


55 . 6 


Use f ul 


71 




80 


67 


11. 1 


Inadequate 


29 




20 


33 


27.3 
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Each statistical regression equation gave much 
better performance than persistence. Persistence gave more 
useful information in August, and likewise the regression 
equations gave slightly more useful information in August, 
The f oil owing is an average of the distribution of forecasts 
for the three-month independent sample by equation category: 





Jul-Aug-Sep 


Jul-Aug-Sep 


Indiv . 


Indiv . 


Equa t ion : 


Recurving 


East-West 


Monthly 

E-W 


Monthly 
Re c . 


Accurate 


47 . 7 


48.3 


42 . 6 


44 . 0 


Adequate 


77.0 


74.7 


72 . 3 


70 . 3 


Useful 


87 . 7 


87 . 7 


86 . 7 


84 . 3 


Inadequate 


12 . 3 


12.3 


13.3 


15 . 7 



In general, the regression equations derived 
from the combined three-month dependent sample gave better 
results than the equations derived from individual months, 
suggesting that equations derived from a larger sample are 
more accurate. In addition, the combined data equation is 
easier to apply. There V7as little difference in results 
whether using an equation derived from east-west or recurv- 
ing data, indicating either equation could be applied with 
equal success. As a matter of convenience, the equation 
from combined monthly data should be used, thus eliminating 
the need to use a different equation for each month and for 
storms which extend into the next month. 

b. Five-predictor equations derived from data with- 
in latitude bands 

The data were next evaluated using the regres- 
sion equations developed from storms in latitude bands. 
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TABLE VIII 



Acceptability distribution of 24-hour intensity forecasts 
in percent of sample size for independent data in latitude 
band 10-19. 9N. 



a. July Data - Sample Size 69 



Equation : 


July 


Jul-Aug-Sep 


Persistence 




10-19 . 9N 


10-19 . 9N 




Accurate 


48 


55 


23 


Ade qu at e 


81 


83 


51 


Useful 


93 


93 


85 . 5 


Inadequate 


7 


7 


14 . 5 


t 


b . 


August Data - Sample Size 75 


Equation : 


Augus t 


Jul-Au g-Sep 


Persistence 




10-19 . 9N 


10-19 . 9N 




Accurate 


27 


47 


36 


Adequate 


58 


76 


56 


Useful 


75 


87 


81 


Inadequate 


25 


13 


19 




c . 


September Data - Sample 


Size 140 


Equation: 


September 


Jul-Aug-Sep 


P er si s t en ce 




10-19 . 9N 


10-19 . 9N 




Accurate 


37 


37 


32 


Adequate 


65 


69 


54 


Useful 


84 


86 


72 


Inadequate 


16 


14 


28 




d'. 


Average for three-month 


period 


Equation : 


Individual Jul-Aug-Sep 


Persistence 




Months 


10-19 . 9N 




Accurate 


37 . 3 


46 . 2 


- 30.3 


Adequate 


68 


76.0 


53 . 7 


Useful 


84 


88.7 


79.5 


Inadequate 


16 


11 . 3 


20 . 5 
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TABLE IX 



Acceptability distribution of 24-hour intensity forecasts 
in percent of sample size for independent data in latitude 
band 20-29,9N, 

a, July Data - Sample Size 92 

Equation: July Jul-Aug-Sep Persistence 

20-29. 9N 20-29. 9N 



Accurate 


39 


40 


31.5 


Adequate 


70 


70 


55.5 


Useful 


85 


85 


68.5 


Inadequate 


15 


15 


31.5 




b. August 


Data - Sample Size 70 


Equation : 


August 
20-29 . 9N 


Jul-Aug-Sep 
20-29 . 9N 


Persistence 


Accurat e 


73 


53 


40 


Adequate 


90 


87 


74 


Useful 


97 


94 


78.5 


Inadequat e 


3 


6 


21 . 5 




c. September Data - Sample 


Size 112 


Equation: 


September 
20-29 . 9N 


J ul-Aug- Sep 
20-29 . 9N 


Persistence 


Accurate 


47 


50 


25 


Adequ a t e 


74 


78 


46 . 5 


Useful 


87 


93 


61 


Inadequate 


13 


7 


39 




d. Average 


for three month 


period 


Equation : 


Individual 

Months 


Jul-Aug-Sep 
20-29 . 9N 


Persistence 


Accurate 


53 


47.7 


32 . 2 


Adequate 


78 


78.3 


58 . 7 


Useful 


89 . 7 


90 . 7 


69 . 3 


Inadequate 


10.3 


9.3 


30.7 
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Due to the lack of sufficient independent data in the 
bands 0-9.9N and 30-39o9N, the equations in these latitude 
bands were not evaluated. The approach used to evaluate 
each equation V7as to apply individual monthly equations and 
also combined July-Augus t-S ep tember equations to monthly 
data in a given latitude band. Tables VIII and IX are ac“ 
ceptability distributions by months for 24-hour intensity 
forecasts in latitude bands using statistical regression 
equations from Tables XII and XIII, in Appendix A and also 
using persistence. Each equation performed better than 
persistence in 24-hour intensity forecasts except the equa- 
tion derived from August data within 10-19. 9N. Regression 
equations to forecast 24-hour intensity for July storms 
gave 93 percent useful information within 10-19. 9N and 85 
percent useful information within 20-29. 9N (Table VIII). 
However, in August and September storms, regression equa- 
tions gave more useful information vrithin 20-29.9N than in 
the 10-19. 9N band. Forecasts using regression equations 
gave more overall useful information within the 20-29.9N 
band, as seen by comparison of part (d) in Tables VIII and 
IX. Considering only persistence, more useful information 
was provided within the 10-19. 9N band. Persistence fore- 
casts yielded more useful information in July and August 
than in September, regardless of latitude. 

Comparison of the distribution of useful infor- 
mation given in Tables VII, VIII and IX with the correla- 
tion coefficient (R) in Table VI provided some very 
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interesting results. The correlation coefficient was 0.74 
for the Jul~Aug"Sep east-west and recurving equations in 
Table VI, part (a). Indeed the amount of useful information 
in Table VIII, part (a) was almost equal, being 86 and 84 
percent respectively for the two equations. From Table VI, R 
for the July east-west equation was 0.73 and for the July re 

curving equation 0.53. As expected, the useful information 

decreased to 83 and 76 percent, respectively, for the two 
$ 

equations. This general trend was consistent throughout 
most of the comparisons. However, the actual magnitude of 
R when considering persistence compared to equations did 
not appear to be a good measure of usefulness. The correla- 
tion coefficient (R) of persistence for August data within 
10-19 . 9N was 0.10 (from Table VI, part (b)), compared to R 
for the Jul-Aug-Sep 10-19. 9N equation of 0.62. Considering 
the distribution of useful information in Table VIII, part 
(b), persistence yielded 8]. percent and the equation yielded 
87 percent. Although there was a marked difference in R, 
the amount of useful information was relatively close. The 
evaluation scheme examining R considers only absolute 
error and applies equal weight to an error regardless of 
the intensity of the storm, and focuses on normalized 
intensity changes, whereas the JTWC method focuses on 
actual intensity. This may account for the differences in 
the distribution of useful information in the FWC/JTWC 
method compared to the correlation coefficients. 
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Considering both latitude bands, 10-19. 9N and 
20-29. 9N, the combined equation derived from July, August 
and September data produced equal or more useful informa- 
tion than equations derived from a single month. The fol- 
lowing is a combined average of the distribution of the 



equations for the two latitude bands 10-19. 9N and 20-29N: 



Equa t i on 


Jul-Aug-Sep 


Individual 

Monthly 


Persistence 


Accurate 


46.9 


45 . 2 


31.2 


Adequate 


77 . 2 


73.0 


56 . 2 


Useful 


89.7 


86.8 


74.4 


Inadequate 


10.3 


13.2 


25 . 6 



The overall averages indicated that two equa- 
tions, namely Jul-Aug-Sep 10-19. 9N and 20-29. 9N, could be 
used for predicting future 24-hour cliange in intensity for 
any storm in the respective latitude band, regardless of 
the month of occurrence. These equations yield significant- 
ly better results than linear extrapolation and slightly 
better results than the individual monthly equations. The 
combined equations from Jul-Aug-Sep data in the latitude 
bands produced 89.7 percent useful information compared to 
87.7 percent useful information from the combined data 
equations in east-west or recurving data. Thus only two 
five-parameter equations are required for forecasting 24-hour 
intensity of tropical cyclones in the months of July, August 
or September. One equation will apply to storms within 
10-19. 9N and the other for storms within 20-29. 9N. Each 
equation gives significantly more useful data than 
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persis teiice , and thus can be expected to perforin better than 
the present methods now being used by FWC/JTWC- As a storm 
moves put of these two latitude bands or occurs in other 
months, a variety of equations are available in Appendix A 
which will yield better results than simple extrapolation, 
c. Ten-predictor equations 

Equations using ten predictors were derived 
from the same data categories as the f ive -pr edi c t or equa- 
tions. Each ten-predictor equation contained the five-para- 
meters which made up the corresponding five-predictor 
equation as well as five additional parameters. Most of the 
five additional parameters were of statistical significance 
at the 5% level when considering the F-value upon entry. 
Although some ten-predictor equations contained predictors 
which were not significant at the 5% level, the overall 
F-value of each ten-predictor equation remained above the 
critical F-value for the 5% level. 

Each ten-predictor equation provided more useful 
information than persistence when forecasting 24-hour in- 
tensity, However, considering the distribution of accept- 
able information from ten vice five-predictor equations, 
there was little difference. A comparison of the percent 
distribution of acceptable information from Jul-Aug-Sep 
combined data ten-predictor equations within 10-19. 9N and 
20-29. 9N verses the corresponding five-predictor equation 
follows : 
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Independent 

Data 



Ten-Predictor 

Equation 

Acc . Ade . Use . 



Five-Pr ed ic t o r 
Equat ion 

Acc. Ad e . Use 



Jul 10-19. 9N 
Aug 10-19. 9N 
Sep 10-19. 9N 



60 

46 

44 



83 

72 

72 



93 

88 

82 



55 

47 

37 



83 

76 

69 



93 

87 

86 



Jul 20-29. 9N 
Aug 20-29. 9N 
Sep 20-29. 9N 



43 
63 

44 



73 

91 

73 



87 

96 

94 



40 

53 

50 



70 

87 

78 



85 

94 

93 



As seen from the above, there was very little 



to be gained by using a ten-predictor equation vice a five- 
predictor equation, and certainly a shorter equation is much 
more desirable from an application point of view. Ten- 
predictor equations to predict actual 24-hour future inten- 
sity are presented in Appendix B with an evaluation of 
acceptable information in Table XVI. 
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IV. SUMMARY AND CONCLUSIONS 



Statistical regression equations to predict future 24- 
hour change in intensity of tropical cyclones were derived 
from a history file of tropical cyclone data stratified by 
monthly periods into the following categories: (a) East- 

west moving storms, (b) Recurving storms, (c) All storms 
within latitude bands 0-9.9N, 10-19. 9N, 20-29. 9N and 30- 
39. 9N. Storms were eliminated which (1) had an initial 
northeasterly heading which continued to move to the north- 
east, (2) had a genesis west of 125E, (3) had a first 

observation reporting a maximum intensity greater than 65 
knots, (4) looped, i.e. formed a complete loop in the track. 
All data were considered over open ocean only. The depen- 
dent sample covered a ten-year period (1960-1969), and the 
independent test sample covered a five-year i^eriod (1955- 
1959) . 

Analysis of forecasts using the independent data sample 
indicated that the statistical regression equations could 
be used to determine future 24-hour changes in intensity of 
tropical cyclones with forecasts which are significantly 
more accurate than current methods based mainly on persis- 
t ence . 

Although typhoon intensity, size, and speed of movement 
varies with month (Liechty 1972), only two five-predictor 
equations were required to forecast 24- hour intensity 
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changes in the majority of July through September storms in 
the western North Pacific. The two equations predict inten- 
sity changes of tropical storms and typhoons occurring from 
July through September within a given latitude band. Since 
none of the past 24-hour predictors were selected for these 
equations, (Appendix A, Tables XII (g) and XIII (f)) only 
12 hours of history are required before a forecast can be 
made. The five predictors are readily available to the fore- 
caster, and v/ithin minutes, a forecast can be made using the 
equations. If tropical cyclones or typhoons occur outside 
the area or periods covered by these two equations, a variety 
of additional prediction equations from which the forecaster 
may choose are included in Appendix A. 

The statistical regression equation to predict future 
24-hour intensities of typhoons must be considered as only 
one additional tool to aid the forecaster. It must be 
remembered that the values of parameters required by these 
equations may contain errors, thus somewhat degrading the 
accuracy of the equations. The only real test will be with 
operational data rather than post-season smoothed, best 
track data. However, these equations, when combined with 
experience, synoptic data, subjective modifications, and 
climatological guides, should greatly improve the accuracy 
of present methods of typhoon intensity forecasting. 
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APPENDIX A 



Statistical regression equations using five parameters 
to predict future 24-hour change in intensity in knots per 
24 hours are presented. Tables X and XI are the monthly 
equations derived from east-V7est and recurving data respec- 
tively. Tables XII and XIII are the monthly equations de- 
rived from data within latitude bands 10-19. 9N and 20-29. 9N 
respectively. Table XIV presents the equations derived 



from data v/ithin 0-9 . 9N and 30-39. 9N. 



TABLE X 



Regression equations derived from east-west data to predict 
future 24-hour change in intensity (kt/24hr). 



(a ) J uly 
Predictor Z 
CONSTANT 
1 

10 

50 

52 

54 

(c) September 
Predictor Z 
CONSTANT 
2 

11 

17 

39 

49 

(e) November 
Predictor Z 
CONSTANT 
1 
14 
17 
19 

50 



Coefficient 
42 . 25240 
-0.99618 
-0.65633 
-0.00996 
-46 .59335 
0.85886 



Coefficient 
-88 . 32253 
0.41633 
0.58877 
0.72737 
1.69331 
-0.00402 



Coefficient 
1024.77002 
-1.45327 
-3 . 11491 
1.14865 
-10.17075 
--0 . 00869 



(b) August 
Predictor Z 
CONSTANT 
2 

10 

13 

50 

55 

(d) October 
Predictor Z 
CONSTANT 
13 
17 
42 
46 
50 

(f) December 
Predictor Z 
CONSTANT 
10 
22 
44 
46 
49 



Coefficient 
-60 . 53868 
0.86914 
-0 .68542 
-1.12739 
-0 . 00928 
2 . 46833 



Coefficient 
-23.36911 
2 . 06295 
1.43920 
0 .81278 
-1.10552 
- 0 . 02001 

through June 
Coefficient 
7 .25270 
-1 . 15147 
-3 . 81075 
0 . 39751 
-0.23558 
-0 . 00099 



(g) Jul-Aug-Sep 
Predictor Z Coefficient 



CONSTANT 

1 

2 

10 

17 

50 



-22 . 35187 
-0 . 87414 
0.36726 
-0.56878 
0.48592 
-0.00566 
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TABLE XI 



Regression equations derived from recurving data to predict 
future 24-hour change in intensity (kt/24hr). 



(a) July 




(b) August 




Predictor Z 


Coefficient 


Predictor Z 


Coefficient 


CONSTANT 


456 .68408 


CONSTANT 


18 . 12944 


10 


-0 .85937 


1 


-1 . 31838 


15 


-0.14767 


10 


-0.59213 


39 


1.11113 


11 


0. 6185.8 


53 


-0.66595 


17 


0 . 32553 


54 


1.20949 


49 


-0.00481 


(c) September 




(d) October 




Predictor Z 


Coef f icient 


Predictor Z 


Coef f icient 


CONSTANT 


48 . 62646 


CONSTANT 


48 . 51552 


1 


-1 . 17217 


1 


-1 .26156 


10 


-0.22306 


10 


-0 . 68253 


17 


0.40493 


17 


0 .40350 


19 


-4 . 64079 


49 


-0.00197 


51 


188 . 81778 


52 


-77 . 76389 


(e) November 




(f) December 


through June 


Predictor Z 


Coefficient 


Predictor Z 


Coefficient 


CONSTANT 


23.45758 


CONSTANT 


-3 .93260 


1 


-1 .20643 


1 


-2 . 30894 


5 


0.03156 


10 


-1 . 16680 


10 


-0 . 85769 


11 


0 .39558 


42 


-1.59244 


13 


1 . 30513 


50 


-0 . 00305 


49 


-0 . 00369 


(g) Jul-Aug-Sep 






P red i c t o r Z 


Coef f ic ien t 






CONSTANT 


43.25537 






1 


-0.99846 






10 


-0.35953 






11 


-0.20295 






17 


0.43457 






19 


-3.10880 
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TABLE XII 



Regression equations derived 


from data within 10-19. 9N to 


predict future 


24-hour change 


in intensity 


(kt/24hr) . 


(a) July 




(b) August 




Predictor Z 


Coefficient 


Predictor Z 


Coefficient 


CONSTANT 


381 .99585 


CONSTANT 


-862 . 32739 


1 


-0. 83874 


10 


-0 .63809 


14 


1.35827 


14 


2 .52380 


18 


0 . 38930 


17 


0.88987 


50 


-0 . 01267 


49 


-0 . 01949 


52 


-64.10526 


51 


2429 . 88159 


(c) September 




(d) October 




Predictor Z 


Coef f icient 


Predictor Z 


Coe f f icient 


constant' 


2 .19348 


CONSTANT 


-0.28871 


10 


-0.40875 


8 


3.20133 


44 


-0.25284 


10 


-0 .76968 


46 


0 . 86847 


11 


0 .42675 


50 


-0.00940 


17 


0 . 52139 


54 


1.52462 


49 


-0.00426 


(e) November 




(f) December 


through June 


Predictor Z 


Coef f icient 


Predictor Z 


Coe f f i c ien t 


CONSTANT 


-31.20357 


CONSTANT 


37.49117 


1 


-2 . 32222 


1 


-2 .25802 


2 


0. 58671 


10 


-1.10579 


10 


-0.39322 


22 


-1 . 63809 


42 


0 .29826 


45 


-0.15726 


49 


-0.00134 


49 


-0 . 00121 


(g) Jul-Aug-Sep 






Predictor Z 


Coefficient 






CONSTANT 


-44.62602 






2 


0.38688 






10 


-0.48747 






17 


0 .58992 






46 


0.63176 






50 


-0.00867 
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TABLE XIII 



Regression equations derived from data within 20-29. 9N to 
predict future 24-hour change in intensity (kt/24hr). 



(a) July 
Predictor Z 
CONSTANT 
10 
16 
21 

50 

54 

(c) September 
Predictor Z 
CONSTANT 
1 
7 

10 

51 

55 

(e) November 
Predictor Z 
CONSTANT 
5 

10 

18 

42 

50 



Coefficient 
372.08813 
-0.76135 
-1.16243 
-4.78572 
-0 . 00632 
0 . 94574 



Coef f icient 
16 . 51608 
-0 .63702 
1.58331 
-0.27586 
-207 . 23415 
5.46685 

through June 
Coefficient 
2.23611 
0 . 03680 
-0.20631 
0.24642 
-2. 66799 
-0.00538 



(b) August 
Predictor Z 
CONSTANT 
1 
2 

10 

19 

49 

(d) October 
Predictor Z 
CONSTANT 
1 
2 
3 

11 

49 



Coef f icient 
39.44713 
-2.04376 
0 .24200 
-0.59548 
-7 . 04820 
-0 . 00166 



Coefficient 
47 . 95605 
-5 .54319 
0 . 40305 
0 . 03157 
0.57115 
-0 . 00417 



(f) Jul-Aug-Sep 



Predictor Z 
CONSTANT 
10 
19 
49 
40 
55 



Coefficient 
6 .56849 
-0 . 45188 
-3 . 30801 
-0 . 00063 
-0 . 00373 
4 .98910 
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TABLE XIV 



Regression equations derived from data within 
30-39. 9N to predict future 24-hour change in 
(kt/24hr) . 



a. Jan through Dec within 
Predictor Z Coefficient 



CONSTANT 

8 

10 

11 

46 

49 



-2 . 86331 
-3.37220 
-1.45857 
0.88116 
-1.87045 
-0 . 00389 



b. Jan through Dec within 
Predictor Z Coefficient 



CONSTANT 

15 

16 
19 
22 
50 



-172 .37785 
-0.28473 
0.67054 
-2.12016 
0.48649 
-0.00462 



0-9. 9N and 
intensity 

0-9. 9N 
30-39 . 9N 
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APPENDIX B 



Ten-predictor regression equations were derived to 
predict actual future 24-hour intensity of tropical cy- 
clones using the data in Section II and the methods of 
Section III, Two equations were derived from combined 
July, August and September data within 10-19, 9N and 20- 
29, 9N, ahd are presented in Table XV, These ten-predic- 
tor equations were evaluated v/ith the forecast verifica- 
tion scheme employed by FWC/JTWC> and the results are 
found in Table XVI, The distribution of useful informa- 
tion obtained from these equations which predicted 
actual future 24-hour intensity is very similar to the 
distribution obtained from forecasts using ten-predictor 
equations v;hich predicted future 24-hour changes in 
ill tens i ty . 
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TABLE XV 



Ten-predictor regression equations derived from combined 
July, August, and September data to predict actual future 
24-hour intensity (in knots) of tropical cyclones within 
latitude bands . 



(a) 10-19. 9N 
Predictor Z 
CONSTANT 

2 

3 

10 

11 

15 
17 
19 
39 

49 

50 

(b) 20-29. 9N 
Predictor Z 
CONSTANT 

1 

2 

10 

11 

1 6 
17 
19 

50 

51 
55 



Coef f icient 
-73.9 

0.56236 
0.03201 
-0.42194 
1.50945 
-0.18754 
0.40647 
1.92792 
0 .42303 
-0 . 00303 
-0 . 00391 



Coefficient 
141.1 
-1.25443 
-0 . 25864 
-0.47648 
1 . 89980 
-0.45965 
-0 . 32026 
-2 .91058 
-0 . 00449 
-895, 2 

6.91715 



t 
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TABLE XVI 



Acceptability distribution of information in percent of 
sample size from 24-hour future intensity forecasts 
within latitude bands using ten-predictor equations. 

(a) Jul-Aug-Sep 10-19. 9N equation to predict 
actual intensity 



Data v/ithin . 


•k 




10-19 .9N 


Jul (69) 


Aug (75) 


Sep 


Ac cur a t e 


61 


46 


46 


Adequa t e 


81 


73 


74 


Useful 


93 


88 


82 


Inadequate 


7 


12 


18 


Jul-Aug-Sep 10-19. 9N equation 
intensity change 


to predict 




Data within ^ 


k 




10-19 . 9N 


Jul (69) 


Aug (75) 


Sep 


Accurate 


60 


46 


44 


Adequate 


83 


72 


72 


Useful 


93 


88 


82 


Inadequa t e 


7 


12 


18 


Jul-Aug-Sep 20-29. 9N equation 
actual intensity 


to predict 




Data within 






20-29 , 9N 


Jul (92) 


Aug (70) 


Sep 


Accurate 


hh 


67 


4 6 


Adequate 


6 8 


90 


73 


Useful 


82 


97 


92 


Inade qua t e 


18 


3 


8 


Jul-Aug-Sep 20-29. 9N equation 
intensity change 


to predict 




Data within . 






20-29. 9N 


Jul (92) 


Aug (70) 


Sep 


Accurate 


43 


63 


44 


Ade qu ate 


73 


91 


73 


Useful 


87 


96 


94 


Inadequate 


13 


4 


6 



^Denotes sample size 
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